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ABSTRACT 

A model is presented in which both Fanaroff and Riley class I and II extragalactic jets are produced 
by magnetized accretion disk coronae in the ergospheres of rotating black holes. It employs a hybrid 
version of the Blandford-Payne and Blandford-Znajek magnetohydrodynamic mechanisms (similar to 
the Punsly-Coroniti model, with the addition of a metric-shear-driven dynamo) and a generalized form 
of the magnetic switch, which is shown to be the MHD analog of the Eddington luminosity. While the 
jets are produced in the ergospheric accretion disk itself, the output power still is an increasing function 
of the black hole angular momentum. For high enough spin, the black hole triggers the magnetic switch, 
producing highly-relativistic, kinetic-energy-dominated jets instead of magnetic-energy-dominated ones 
for lower spin. The coronal mass densities needed to trigger the switch at the observed FR break power 
are quite small (~ 10~ 15 gcm~ 3 ), implying that the source of the jet material may be either a pair 
plasma or very tenuous electron-proton corona, not the main accretion disk itself. 

The model explains the differences in morphology and Mach number between FR I and II sources and 
the observed trend for massive galaxies (which contain more massive black holes) to undergo the FR I/II 
transition at higher radio power. It also is consistent with the energy content of extended radio lobes 
and explains why, because of black hole spindown, the space density of FR II sources should evolve more 
rapidly than that of FR I sources. 

A specific observational test is proposed to distinguish between models like this one, in which the FR 
I/II division arises from processes near the black hole, and models like Bicknell's, in which the difference 
is produced by processes in the host galaxy's interstellar medium. If the present model is correct, then 
the ensemble average speed of parsec-scale jets in sources distinguished by their FR I morphology (not 
luminosity) should be distinctly slower than that for sources with FR II morphology. The model also 
suggests the existence of a population of high-redshift, sub-mJy FR I and II radio sources associated 
with spiral or pre-spiral galaxies that flared once when their black holes were formed but were never 
again re-kindled by mergers. 

Subject headings: black hole physics — galaxies: jets — galaxies: nuclei — hydrodynamics — MHD — 
quasars: general — radio continuum: galaxies — relativity 



1. INTRODUCTION 

Recently many authors have pointed out that the speeds 
of jets associated with different types of compact object 
(pre-main-sequence star, white dwarf, neutron star, black 



Radio jets on the kiloparsec scale also differ dramatically 
in speed, and these differences apparently translate into a 
morphological difference. Subsonic or transonic jets with 
velocities less than ~ 0.6c are believed to be respons ible 



hole) arc approximately equal to the Kcplcrian or escape 
apood from the ourfacoo of thooo compact objocto. — How 



for producing the Fanaroff and R i ley Class I sou rces ( Fa- 



naroff fc Riley 1974| ; [Bickncll 198g [Bickncll 1995|) . The FR 



ever, some caution must be exercised in taking this equiv- 
alence too literally. A good example of this is the extra- 
galactic radio source class. All are believed to be pow- 
ered by accretion onto black holes, from whose surface 
(the horizon) the escape velocity is exactly c, the speed of 
light. However, while all extragalactic jets appear to be 
at least mildly relativistic, the speeds obtained from Very 
Long Baseline Interferometry (VLBI) observations range 
from (3j = Vj/c ~ 0.1 in Centaurus A ( Lingay et al. 1998 ) 
to jj = (1 - (3f)- 1/2 ~ 10 in 3C 345 ( |Vcrmculcn 1996| ) 



In terms of the proper velocity Uj — Jj/3j, this represents 
a range of two orders of magnitude in jet speed. Clearly, 
even before traversing the main part of the interstellar 
medium, jet speeds are strongly affected by the details of 
the acceleration and collimation process near or in the ac- 
creting black hole system. 



I sources are distinguished by generally low radio power, 
often complex morphology (S-shaped, C-shaped, etc.), and 
diffuse lobes that have their bright region close to the as- 
sociated galaxy — all consistent with transonic flow that 
decelerates into subsonic flow. Class II sources generally 
are much more powerful, have only a straight morphol- 
ogy, and have their brightest point at the far ends of the 
dumbbell-shaped source. They are thought to be produced 
by highly supersonic jets that terminate in a strong shock 
where the high-speed flow strikes the interstellar medium. 
A mechanism for determining the speed of jets at their 



(1997a) 



site of acceleration was suggested by Meier et al 
hereinafter MEGPL. Termed the "magnetic switch" be- 
cause of the two distinct slow and fast jet states, this 
mechanism operates in the magnetized corona of the black 
hole accretion disk itself. Fast (highly relativistic) jets are 
produced when the coronal plasma is unbound by the mag- 
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netic field, and slow (trans-relativistic) jets at roughly the 
disk escape velocity are produced when the corona is ini- 
tially bound. This simple physical process is a good candi- 
date for explaining the difference between the two Fanaroff 
& Riley classes, because it explains how jets with radically 
different speeds could occur in otherwise similar sources. 
However, as discussed in section || below, the magnetic 
switch, triggered solely by the magnetic field strength (or 
even by coronal density), fails to explain why these mor- 
phological differences have rather long lifetimes (10 6-7 yr). 

One way out of this dilemma is to assume that all radio 
jets for a given galaxy mass indeed do start out in the nu- 
cleus with about the same average speed, but, depending 
on local galaxy properties (which do have long lifetimes), 
are affected differe ntly by deceler ation processes in the in- 
terstellar medium ( Bickncll 1995| ). Indeed, jet deceleration 
has been inferred from variations in the jet-to-c ounterjet 
ratio in observations of FR I sources flLaing 1996| ). In this 
paper, however, we propose an alternate method for main- 
taining long-lived FR I or FR II morphology — one that 
preserves the magnetic switch model, along with all of its 
attractive properties, and suggests a reason why some jets 
may be more prone to deceleration than others. Since the 
power in these MHD outflows also depends on the angular 
velocity of the magnetic field, it is suggested that the dif- 
ference in radio power among galaxies of the same mass is 
due to different magnetic field rotation rates in their cen- 
tral engine, caused by different spin rates of their central 
black holes. Th is has been suggested before ( |Baum, Zirbcl " 



& O'Dea 1995) based on observations alone, but this paper 
derives such a behavior directly from a generalized theory 
of the magnetic switch. As we show below, black holes of 
different spin should produce jets of different power and 
velocity, and a sharp FR I/II-like transition is expected to 
occur. Unlike the ISM model, this model predicts a dif- 
ferent mean jet velocity in the central parsec for the FR 
I and FR II classes, yielding a definite observational test 
for distinguishing between the two. 

Note that the magnetically-switched, rotating black hole 
model depends critically on the jet luminosity being a di- 
rect function of the black hole angular momentum — in 
direct conflict with the conclusion s of Livio, Ogilvie, & 
Pring 



(1999) 



Therefore, section |3.2| spends some time 
on this issue and shows that, while the basic tenets of 
these authors are correct, nevertheless, if the accretion disk 
within the ergosphere is taken into account, a Blandford- 
Znajek-like (BZ) expression for the jet power is recovered 
- even though the jet is ejected from the accretion disk 
itself by a Blandford-Payne-likc (BP) mechanism. The jet 
production model is similar t o the Punsly-Coroniti (PC ) 
ergospheric wind mechanism ( Punsly fc Coroniti 1990b| ). 
Furthermore, due to an expected metric- shear- driven dy- 
namo (see the Appendix), the poloidal magnetic field in 
the ergospher e is expected to have a larger valu e than that 
suggested by Livio, Ogilvie, fc Pringle (1999)| , leading to 



radio jet powers comparable to those observed. 

The outline of the paper is as follows. Section [2] briefly 
reviews the classical magnetic switch model and the rea- 
sons for its failure to explain the FR I/II break. Section 
H presents the generalized magnetic switch (including ro- 
tation as a parameter) and its role in the rotating black 
hole model for jet production and for the FR I/II transi- 
tion. Section discusses applications of the model to re- 



cent observations of radio sources, with observational and 
theoretical tests of the model proposed in section [| 

2. THE MAGNETIC SWITCH FOR NORMAL COMPACT 
OBJECTS AND NON-ROTATING BLACK HOLES 

2.1. The Classical Magnetic Switch 

In the original magnetic switch model presented in 
MEGPL, the flow behavior depends on the ratio of 
the Alfven velocity in the inner disk corona Vao = 



,1/2 

a/2 



cO 



to the escape velocity there V e . 



(2GAf/i? ) 1/2 , where B p0 is the strength of the poloidal 
magnetic field protruding from the disk into the corona at 
radius i?o, Pco is the density of the corona there, M is the 
mass of the black hole, and G is the gravitational constant. 
When v = Vao/VcscQ < 1, gravity dominates magnetic 
forces (the plasma is bound initially) and the character 
of the flow is similar to a Parker wind (with collimation) : 
azimuthal magnetic field coiled in the corona uncoils up- 
ward to slowly accelerate the flow through a series of crit- 
ical points. Depending on the parameters, the resulting 
jet speed is approximat ely equal to the escape speed at 
the base of the outflow (Kudoh & Shibata 1996), which is 
somewhat interior to the last stable orbit in the relativis- 



tic Schwarzschild black hole case (Koide et al. 1998a). 
When v > 1, magnetic forces dominate gravity and, for 
sufficiently high v, render the latter negligible. (The coro- 
nal plasma is unbound.) As in pulsar wind theory, the 
outflow velocity then is determined not by gravitational 
effects, but by the magnetization parameter 



*o = Vi Rl/(AV M cR 2 L ) 



(1) 



which is a measure of the amount of work done on the in- 
jected plasma; here, Vio is the velocity with which matter 
is injected, Rl = c/Q is the radius of the light cylinder, 
and fl is the angular velocity of the poloidal magnetic field 
at i?o- For relativistic flows, jj ~ ctq appears to be a g ood 



approximation for the final M HD wind/jet velocity (Ca- 
menzind 1989] ; Henrikscn 1989 ) 



A transition between the two types of outflow then is 
induced by increasing the magnetic field strength B p q un- 
til v > 1, whereupon significantly higher jet speeds occur. 
This magnetic switch has several properties that make it 
an attractive model for the difference in morphology be 



tween class I and cl ass II extragalactic radio sources ( Fa 



naroff fc Riley 1974 ): 1) the transition is sharp, occurring 



with increasing jet power; 2) jets produced when v < 1 are 
transonic, while those for v > 1 are highly supersonic, and 
3) even statistical behavior observed in large samples of 
radio galaxies o f different power an d optical m agn itude is 
reproduced (see Meier et al. 1997b and Section 4.1 below). 



2.2. Failure of a Change in Magnetic Field Strength to 
Explain the FR I/II Break 

However, as presented in MEGPL, the original magnetic 
switch, triggered by either increasing the coronal magnetic 
field or decreasing the coronal density, fails to explain the 
FR I/II break, chiefly because it fails to explain why ra- 
dio sources appear to maintain their FR morphology for 
long periods of time — at least as long as the flow time 
from galactic center to lobe (10 6-7 yr). There is no rea- 
son to expect the strength of the magnetic field threading 
disk and hole (or the coronal density) to remain fixed at 
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high or at low values over such long time scales. Indeed, 
one expects that these quantities will be distributed fairly 
equally about a mean, and on average will be the same 
for galaxies of the same mass. The argument is as fol- 
lows. The structure of an accretion disk is a function of 
the black hole mass rrig = M / 10 9 Mq, the accretion rate 
M, the viscosity parameter a, and position in the disk 
scal ed to the Schwarzschild r adius (Bhakura & Sunyae\ 
1973; Novikov fc Thornc 1973 ). Now, there is growing ev- 
idence that the black hole mass, even in presently inactive 
objects, is approximately proportional to th e mass of the 
centr al spheroid al bulge of stars (see, e 
Richs .one 1995 ) such that 



Kormendy & 



log m 9 



-0.5Mb - 10.5 



(2) 



where Mb is the absolute blue optical magnitude of the 
bulge. In addition, while the instantaneous accretion rate 
will vary dramatically on short time scales (years), be- 
cause the hole is fed by galactic stars and gas, over the 
very long term the average accretion rate also should be 
roughly proportional to the bulge mass (a perhaps neces- 
sary condition for equation to be valid). Therefore, to 
first order, we expect the average accretion rate 

< to > ee < M > I Mem =< M > /(1.4 x 10 26 gs" 1 m 9 ) 

(3) 

to be a constant, proportional to the average amount of 
matter shed per unit time into the nucleus and onto the 
hole by a galactic bulge of a given optical magnitude, 
divided by the mass of the black hole in that nucleus. 
(Note that the definition of m here does not include an 
efficiency of conversion of accreting matter into radiation, 
allowing for maximum values somewhat greater than unity 
[riiEdd = 2.5 — 10], depending on the spin of the black hole. 
It also assumes that the accretion is continuous; to take 
intermittent accretion into account, we will add a duty cy- 
cle parameter below.) Furthermore, d etailed simulations 



of the magneto-rotational instability (Btonc ct al. 1996 



Brandenburg ct al. 1996 ) have shown that a ~ 10 is 



a fairly good approximation for turbulent accretion flows, 
presumably independent of black hole mass to first order. 
And, finally, jets are believed to be ejected from near the 
black hole horizon (i.e., from just inside the last stable or- 



bit, Koide et al. 1998a ), rendering the position parameter 
Ro/Rsch a near-constant value for all objects. Therefore, 
ignoring external influences, over periods of time compa- 
rable to the flow time from nucleus to lobe, all parameters 
that affect the disk structure equations are expected to be, 
at most, functions of the black hole/galactic bulge mass 
only. On average, one expects bulges of the same mass 
to produce the same type of accretion disk in the nucleus, 
with similar magnetic field and coronal density, and, there- 
fore, the same type of jet (either all slow or all fast). It is 
not expected that there will be a spectrum of jet structures 
that maintain themselves for millions of years. Yet the ob- 
servations show that just such a range of radio powers and 
morphologies exists, and that a radio source remains an 
FR I or II for longer than a flow time. 

3. THE MAGNETIC SWITCH IN THE ROTATING BLACK 
HOLE ENVIRONMENT 

In order to lift the degeneracy in expected jet speeds 
and morphologies, we shall take as an additional parame- 



ter the rotation rate of the magnetic field, as determined 
by the black hole rotation rate itself. This section develops 
the model of the accreting rotating black hole system, the 
production of radio jets, and the triggering of the magnetic 
switch. 

In order to compute observable quantities, we must as- 
sume an accretion scenario in the black hole system. For 
definiteness, we choose a rather benign (standard) accre- 
tion disk in which electron scattering dominates free-free 
absorption and gas pressure dominates radiation pressure, 
and we scale the results to a rather low accretion rate 
(to ~ 10~ 3 ). While the details will change depending 
on the exact accretion model and rate, the basic con- 
clusions herein will not, even in the extreme cases of 
near-Eddington accretion (radiation pressure dominated) 
or very low, ad vect ion-dominated accretion flow (ADAF; 
Narayan fc Yi 1995|). 



3.1. The Generalized Magnetic Switch 

The first task is to generalize the magnetic switch theory 
to include rotation as a triggering parameter, in addition 
to magnetic field strength and plasma density. Consider 
a general magnetic rotator of size Rq, poloi dal magnetic 



field Bpp, and angular velocity fi. Following Blandford & 



Payne (1982) , the work done on the plasma trapped in 
the field lines, resulting in an outflowing magnetohydro- 
dynamic (MHD) wind or jet, is 



L 



MHD 



B p0 -Ro^/2 



(4) 



for non-relativistic flow. Now, for a normal Parker-type 
wind, the plasma in the injection region is in hydrostatic 
equilibrium, with only a slight force imbalance outward 
that slowly accelerates the outflow through a critical point 
until it reaches es cape velocity and leaves the system (see, 



Meier 1982[ ). This is also true of MHD winds, al- 



e.g., _ 

though the critical point structure and geometry are more 
complex. However, based on the magnetic switch mech- 
anism described by MEGPL, there should exist a critical 
power, given by the liberation of an escape energy in a 
free-fall time 



Lcrit — E esc0 /Tff — 4:TTp c0 R£ 



( GM 



3/2 



(5) 



such that if Lmhd exceeds L crit , then the plasma becomes 
unbound — even in the injection region itself. That is, the 
outward magnetic and centrifugal forces there strongly ex- 
ceed gravity, accelerating the plasma to greater than the 
escape velocity within Rq and eventually to a much greater 
final velocity The magnetic switch luminosity, then, plays 
the same role in MHD acceleration as the Eddington lu- 
minosity plays in radiative acceleration. For Keplerian 
disks (with Q = 9, K = (GM/i? ) 3/2 ), the condition 
Lmhd > L cr a reduces to the original magnetic switch 
condition Vao > Vesco- However, for a non-Keplerian ro- 
tator of fixed magnetic field strength and plasma density, 
there should exist a critical angular velocity 



a 



V 2 

o escO 
V A0 



(6) 
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such that if f2 exceeds £l C riti then the magnetic switch 
also should be triggered and a fast jet of the type found 
by MEGPL should ensue. 

The author tested this hypothesis by performing numer- 
ical MHD simulations similar to MEGPL, but holding B p q 
fixed and varying the disk rotation speed instead. The 
left panel of Figure [l] shows the resulting jet speed as a 
function of Q. For angular velocities below £l C rit the ter- 
minal jet/wind speed is approximately the escape speed 
at the injection point. When f2 > 0, cr u, the material 
is unbound and quickly accelerated (in a well-collimated 
jet, similar to MEGPL) to a speed such that a large frac- 
tion of the power L mh d is in the form of kinetic energy 
(see the right panel of Figure |l]). It is this strong dif- 
ference in the character of the outflows below and above 
Qicrit t licit is so similar to the FR I/II division: while the 
total power varies smoothly as Lmhd exceeds L cr i t , the 
character of the outflow varies dramatically from a tran- 
sonic, magnetic-energy-dominated jet to a highly super- 
sonic, kinetic-energy-dominated one. 

The details of the behavior of the curves in Figure [l] can 
be understood qualitatively as follows. When £1 < fix , the 
injected plasma is no longer supported fully by rotation. 
Therefore, as soon as it emerges from the disk boundary, it 
dynamically contracts in cylindrical radius, seeking a new 
Keplerian equilibrium. The jet then is ejected from this 
new, smaller, radius R' with a speed similar to the Kep- 
lerian velocity there and a power slightly larger than that 
given by equation (0), because of compressional increase 



in the axial magnetic field component (B' z 



zo 



(xR'n 2 ). 



if n 



is decreased further, the plasma contracts more, increas- 
ing the eventual jet velocity slightly. On the other hand, 
when Q > fix, the reverse occurs, with the plasma ex- 
panding to a new injection radius. As fl increases, the 
value of flcrit calculated from equation (||) changes as well, 
chiefly because only the radial component of the magnetic 
field is important in the expanded state (B' R0 cx R' ~ 
This increases the critical value to fl' crit — fl'x^Lc 



2 IV'ar 



(where V^ R is the radial component of the Alfven veloc- 
ity at the new injection radius), which is greater than the 
value in equation (^|). Finally, well above the magnetic 
switch (when fl >> fl C rit) the well-collimated jet flow 
gives way to a less-collimated radial wind, and the increase 
in outflow speed with rotation speed begins to break down. 
It is not clear whether such a decrease at high fl will occur 
in the relativistic case, where light cylinder effects will be 
important. Relativistic MHD simulations are needed to 
investigate this behavior. 

Some brief comments should be made about boundary 
conditions. In these simulations, the kinematic and mag- 
netic properties of the inflowing matter at the equator (the 
"disk corona" ) are fixed throughout the calculation. How- 
ever, some authors allow these parameters to be deter- 
mined by the solution, particularly the azimuthal and ra- 
dial components of the fi eld at the equator and the inflow 



velocity there (see, e.g., Ustyugova et al. 1999| ). Both of 



these sets of conditions (fixed or self-adjusting) are easy 
to apply numerically, and each is appropriate under dif- 
ferent circumstances. The latter is u seful for generating 



steady state solutions to compare with Blandford & Payne 
(19820's self-similar, semi-analytic models. In this case, 
conditions at the base of the outflow re-adjust themselves, 



subject only to the strength of the vertical magnetic field, 
the density of inflow, and the speed and shear of the disk 
rotation. The system finally reaches a steady state with 
parameters lying on (or near) the Blandford-Payne solu- 
tion surface in (n, A, £o)-space (see BP), with the flow vary- 
ing smoothly from the outer regions all the way back to 
the equatorial boundary. 

On the other hand, fixing the inflow conditions on the 
equator, as we have done here, is perhaps more appro- 
priate for studying the astrophysical situation where the 
magnetic field is anchored in an accretion disk (or other 
type of rotator). In that case the local magnetic field is 
determined by processes internal to the rotator, not by 
the outflow that originates above it. When the boundary 
conditions are fixed, a steady flow still develops above the 
disk, and the azimuthal and radial fields and flow velocity 
in the corona also re-adjust. However, because that result- 
ing outflow does not necessarily match onto the fixed disk 
boundary conditions, there is a transition region between 
the disk and the base of the outflow region, in which the 
magnetic field and velocity dramatically change in magni- 
tude and direction. It is not yet fully clear whether this 
transition region merely provides a way to connect an out- 
flow solution that would have occurred anyway to the disk 
conditions, or whether it plays an active role in determin- 
ing what outflow solution is chosen. What is clear is that 
the values of the parameters chosen on the disk below the 
transition region do matter: the flow changes dramatically 
with the strength of the poloidal and azimuthal field spec- 
ified and, to a lesser extent, with the angle of the poloidal 
field (|Meier, Payne, fc Lind 1996|; MEGPL). Furthermore, 



jet simulations in which a substantial equatorial field (az- 
imuthal plus radial) is specified seem to display more colli- 
mation, presumably because of the strong azimuthal field 
that is initially present or that develops after a few turns 
of the disk. 

3.2. The Rotating Black Hole Model of Jet Production 

The proposed model of jet production from the accret- 
ing, rotating black hole system is ess entially a marriage 
betwee n the disk acceleration model of Blandford fc Payne] 
1982) and the rotating blac k hol e acceleration model of 
Blandford fc Znajck (1977)| and |MacDonald fc Thorne 



1982). The former authors considered the case when 
the magnetic field threads only the Keplerian portion of 
the disk, while the latter examined the situation when it 
threads both hole and disk. In this paper we consider the 
intermediate case when the magnetic field threads both a 
region of the system outside the static limit that does not 
experience appreciable frame dragging and a region closer 
to the black hole that does — the ergosphere. The main 
features of this model are the following: (1) Unlike the 
BZ model, the magnetic field does not have to thread the 
horizon itself in order to extract the black hole rotational 



energy, which avo ids the problems discussed by Punsly & 
Coroniti (1990a), and the primary source of material in 
the jet can be a corona created by the accretion disk it- 
self, not necessarily particles created in spark gaps near 
the horizon. (2) It is unlike the BP model because the 
disk and field there are frame-dragged by the hole rela- 
tive to an observer at infinity, and thus the jet production 
couples well to the black hole spin, extracting rotational 
energy even though it is disk rotation that accelerates the 
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jet. 

This model is very similar to that of PC, who studied 
steady state winds from the ergospheric region of a rotat- 
ing hole. We note here an additional feature of such a 
system that occurs when time dependence is taken into 
account: as shown in the Appendix, when a plasma with 
an appreciable poloidal magnetic field is present in the er- 
gospheric region, the differential dragging of frames will 
increase the azimuthal magnetic field at the expense of 
the hole's rotational energy. If, as a result of dynamo ac- 
tion, this also is accompanied by a corresponding increase 
in the poloidal component, this will have the effect of in- 
creasing the jet power produced near a rotating black hole 
to a value substantially larger than that for a disk around 
a Schwarzschild hole. 

3.2.1. MHD output power 

When a black hole is spinning rapidly, the accretion disk 
can extend well into the ergosphere — the region where 
local observers must rotate with the hole — with equato- 
rial radius R e = 2GM/c 2 . For a nearly-maximal Kerr hole 
(j — » 1), the last stable orbit reaches into this region — 
very close to the horizon at Rh = (1 + yl — j 2 ) GM/c 2 , 
where j = J / (GM 2 / c) is the normalized angular momen- 
tum of the hole (equal to "a/M" in geometric units). Now, 
consider a loop of coronal magnetic field B p q with one foot 
anchored in the ergosphere and one ou tsid e. 1 Recently 



Hayashi, Shibata, fc Matsumoto (1996) and Romanova et 



al. (1998) have shown that even in a Keplerian flow, differ- 
ential rotation in the inner disk can modify closed coronal 
loops, stretching and twisting them into long, virtually- 
open vertical helical structures of the type needed to ac- 
celerate MHD winds and jets. If the differential rotation is 
further augmented by frame dragging, then the work done 
on the field and outflow derives its power not only from the 
Keplerian rotation of the disk, but also from the rotation of 
the ergosphere itself relative to the point where the outer 
foot is anchored (see the Appendix). The system then 
appears to an observer at infinity as a simple magnetic ro- 
tator with size Rq — R e , magnetic field B p q, and angular 
velocity ft > ft e w 0.4ft ff w 0.4j/(2GM/c 3 ), where ft e 
and ft# are the angular velocity of the ergosphere equator 
and hole, respectively. The form of equation (f|) for a rel- 
ativistic rotator 2 yields an MHD power for the accreting 
system of 



Lmhd — B 2 Rq ft 2 j 4c 



1.1 x 10 4S ergs" 



\ 10 5 G 



™2 -2 

m 9 j 



(7) 



(8) 



which is a comparable power to that of the BZ mechanism 
for the same poloidal field strength. 

3.2.2. Disk magnetic field strength 



by equation [Q) and that from the black hole (the BZ pro- 
cess, given by equation^ with ft = ftjy) are, at best, com- 
parable, and, at worst, entirely dominated by the former 
if the hole spin is small. They, therefore, conclude that 
the BZ process will never be important and that the spin 
of the black hole is probably irrelevant for the expected 
electromagnetic power o f the system. This author agrees 



with the basic tenets of Livio, Ogilvie, fc Pringlc (1999) 



that the disk and hole power will be comparable, but not 
with their ultimate conclusion on the unimportance of the 
spin. The key difference is that, in the ergosphere, the 
MHD power of the disk itself (equation ||) should depend 
on the spin of the black hole. In addition, when the hole 
is rotating, a much larger poloidal magnetic field strength 
is allowed than estimated by these authors, yielding black 
hole rotation-driven powers comparable to those observed 
in extragalactic radio sources. 

First, consider the non-rotating Schwarzschild black 
hole case. The disk will be a standard Keplerian one, w ith 
internal disk field given by (Shakura & Sunyaev 1973) 



B d0 

8tt 



apdo 



(9) 



where pdo is the local internal disk pressure. Now the 
source of a coronal loop of field B p0 will be this (largely 
azimuthal) disk field, first protruding into the corona and 
then growing as the loop is sheared and expanded upward. 
However, in a thin Keplerian disk, if the poloidal field is 
allowed to grow to the same order as B^o, then the ex- 
pected MHD power from equation ([|) would be orders of 
magnitude greater than the accretion luminosity 

Lmhd(B p o = B d0 ) _ ac s0 



L„ 



\V R0 \ 

3.8 x 10 3 (a_ 2 ™ 9 ) 1/10 ml 3 1/5 (10) 



where c s o is the sound speed in the disk, Vro is the in- 
ward accretion drift speed, and the notation a_2 repre- 
sents a/10 -2 . (As discussed earlier, we have used a gas- 
pressure, electron-scattering-dominated disk, but the re- 
sults are similar for other types of disks.) There is not 
enough accretion luminosity to support this MHD power 
by factors of several thousand. As B p0 grows to a even 
a small fraction of Bdo, the MHD disk power becomes 
a substantial fraction of L acc . At that point, a signifi- 
cant amount of angular momentum is removed from the 
disk by the outflow, altering the disk structure and en- 
hancing the local accretion rate. In order that the ac- 
cretion draining the inner disk remain in a steady state 
with that feeding it from the outer disk (and not, for ex- 
ample, launch into a series of bursts that each destroy 



Livio, Ogilvie, fc Prmgle (1999) have argued that the or B 



the inner disk entirely), B v p must be self-l imited to about 
the value postulated by Livio, Ogilvie, & 

where 



p0 



Bdo{Ho/Ro) « B, 



do . 



Pringle (1999) 
the half- 



o is 



MHD power from the accretion disk (the BP process, given thickness of the disk at Rq. 

lr This distinction is somewhat arbitrary, as frame dragging takes place outside the ergosphere as well, although at an increasingly smaller 
rate as one moves away from the hole, and it also increases within the ergosphere as one moves toward the horizon. The important point to 
note is that a field loop with its feet at different disk radii will receive a twist similar to that in a normal Keplerian disk, but driven by the 
differential drag 

2 Equation (R 



the mass lost in the jet. 



;ing of frames in addition to any shearing fluid flow within the ergosphere. 
I also can be derived from equation (Q) and the relations "fj = <tq and Lmhd 



"fjMjC 2 , where Mj 



47r R% pco Vio is 
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However, when the central object is a rotating black 
hole, the situation will be quite different. As before, when 
B p o grows to a fraction of Bdo, a significant MHD out- 



flow will develop, exerting a braking torque on the inner 
portion of the disk. Because of the differential frame drag- 
ging, however, even if all the Keplerian angular momentum 
were removed by the outflowing MHD power, the disk ma- 
terial still will be sheared by the metric itself, continuing to 
increase the local magnetic field strength and continuing 
to produce MHD power. The Appendix shows quantita- 
tively how this occurs in the Kerr metric, producing field 
enhancement everywhere in the ergospheric region. Ten- 
sion in this growing, working field then produces a back- 
reaction on the disk plasma falling toward the hole, ac- 
celerating it relative to the rotating frame in a direction 
retrograde to the black hole spin. By the time the magne- 
tized accretion flow enters the black hole, it could have a 
significant negative angular momentum which, when added 
to the hole, would decrease the spin. Therefore, because of 
frame dragging, power flows from hole to disk to field and 
finally to the outflow, exerting a braking torque on the hole 
itself. No longer self-regulated by the accretion flow, B p q 
in the ergospheric region above the disk now can continue 
to grow at the expense of the hole rotational energy to 
order Bdo , given by the d isk solutions with a Kerr metric 
(Novikov & Thornc 1973). For the gas pressure/electron- 
scattcring accretion model assumed here, this value is 

n> id o a , m4fi 1/20 -9/20 . 2/5 /-, , n 

B p q « B d() = 3.4x10 Gq. 2 m Q 1 mj 3 (11) 



Equation (0) then becomes 



L 



MHD 



= 1.3 x 10 47 ergs 



1/10 11/10 . 4/5 -2 



rn 



in 



-3 3 (12) 



dcrski & Sikora 1996 
and" a 



(This differs from equation (16) of [Mo 
by a factor of a, from equation (^), and a factor of ~ 0A Z 
from the slightly smaller ergospheric rotation rate.) While 
Lmhd appears to scale with mg and m in a manner simi- 
lar to an accretion luminosity, it is important to note that 
it does so only because of the scaling on Bdo- The en- 
ergy liberated is derived mainly from the rotation of the 
hole itself and potentially much larger than the accretion 
luminosity (L acc w 10 43 ergs _1 mg rh_s). 

3.2.3. Observed radio power 

To turn Lmhd into an observed radio power, we note 
that, over long periods of time, the central engine may 
generate a jet only a fraction £ of the time. Further- 
more, only a fraction e of this average power is radiated at 
GHz frequencies in a 1 GHz bandwidth and observed with 
standard radio telescopes. From millimeter and optical 
observations of re peated outbu rsts of radio galaxies and 
quasars (see, e.g., Robson 1992), we estimate typical duty 
cycle values to be £ ~ 0.01 — 0.1. The radio efficiency is 
also very uncertain, but we note that for the microquasars 
GRO J1655-40 and GRS 1915+105, estimates of e are of 
order 10~ 5 or so, just for the instantaneous moving jet 
( Meier 1996| ). A canonical value of e ~ 10~ 3 was chosen 
for these calculations by comparing equation ( |l2| ) with the 
maximum observed radio source power (lO^WHz- 1 ) for 

3 Another possible scaling would be with both pdo and the internal 
slope (-0.55) in equation (flSJ) below, but still within the observational 



a radio galaxy that should have a 10 9 Mq hole in its nu- 
cleus (Mb ~ —21.0). This value for e is consisten t with 
total conversion rates of ~ 0.01 — 0.1 (see, e.g., Leahy 
199l[ ), with only a fraction of that radiated in a 1 GHz 



bandwidth. Scaling to these values, the calculated radio 
power for different black hole spins is 



Prad = K Lmhd / 10 9 'Hz 
= 1.3 x lO^WHz" 1 m 



11/10-2 

9 J 



(e_ 3 C-i a 



1/10 - 4/5 
2 m -3 



(13) 



(14) 



where m is still the average outburst accretion rate only, 
not that averaged over the entire duty cycle. 

3.3. The Magnetic Switch as the FR I/II Break 

The final step in the model is to identify the magnetic 
switch as the cause of the FR I/II break. That is, one 
should observe a distinct change in the character of the jet 
produced when Lmhd = L cr n, or when the radio power 
attains the critical value 



P. 



FR 



rad,mod 



eCi c „t/10 9 Hz 



(15) 



To compute an observed break power we need an expres- 
sion for the corona density p C Q at Rg. Of course, coronal 
physics is poorly understood, even for the sun. Neverthe- 
less, we can obtain a rough estimate by assuming that the 
corona is produced locally by the accretion disk, so the 
disk and corona densities should scale together as 3 



PcO 



VPdO 



(16) 



where pdo is the internal disk density at Rq and r\ << 1 
is a parameter which, while not necessarily a true con- 
stant, is assumed to vary sufficiently slowly with the other 
parameters that it does not affect significantly the scal- 
ing of our results with mg, j, and m. The parameter 
rj is expected to be quite small; in the solar corona, for 
example, 770 = P c q/pq{t = 1) ~ 10~ 8 — 10~ n in the 
region between the photosphere and 2i?Q. If the corona 
is a pair plasma, 77 (which is a mass ratio) could be very 
small indeed. Combining equations (||), (|l5|), and (16), 
and again using the gas pressure/electron-scattering cu 
one obtains 



is 



Pfadanod^ 4.8 xl0 25 W Hz" 1 „, 
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• (6_ 3 C-l"-2 7/10 ^ 2/ 3Vll) (17) 



This is the magnetic switch model prediction for the locus 
of the FR I/II break as a function of black hole mass. 

4. APPLICATIONS AND DISCUSSION 

This section applies the model to the FR I/II break 
in the radio-optical plane, shows that the rotationally- 
triggered magnetic switch is consistent with long source 
morphology lifetimes, and addresses the more general 
questions of evolution and radio loud vs. radio quiet 
sources. The important features of the model and its ap- 
plications are shown in figure |[ 

disk temperature T^o (i.e., with Pdo)- This yields a slightly weaker Mr 
errors. 
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4.1. The "OLWL" Relation in the Radio-Optical Plane 



Owen fc Lamg fiy89 T, pwen fc White (1991)| , and 

" " and FR 

type radio galaxies occupy different regions of the 



have shown that FR I 

IT - 

log 

Prad.obs - M r ^a.s plane, where logP radi0fcs is the ob- 
served radio power at 1400 MHz (in WHz -1 ) and 
is the galaxy optical red magnitude measured to a surface 
brightness of 24.5 magnitudes per square arcsecond. In 
particular, the tendency for a galaxy of a given magnitude 
to transition from an FR I to an FR II appears to be a 
strong function of optical magnitude. The dividing line is 
given approximately by 



logP, 



FR 

rad,i 



= -0.66M fl ,24.5 + 10.35 



(18) 



or PrJd,obs x Llp t 5 - The error in the coefficient of M flj2 4.5 
is at least ±0.1 and in the intercept at least ±0.3 — 0.5. 
The Owen-Laing-White-Ledlow (OLWL) relation is a more 
precise statement of the Fanaroff & Riley class division, so 
care must be taken when classifying radio galaxies based 
on their radio power only. There exist FR II radio galaxies 
with quite small powers (occurring in less bright elliptical 
galaxies), as well as FR I objects with very high radio 
power (occurring in very bright galaxies). 

From the magnetic switch model for the FR I/II break, 
it is possible to calculate a theoretical relation for the ob- 
served break power to compare with equation ( fL8| ) . Using 
equation ( p7| ) to define the break, equation (0) to convert 
black hole mass to galaxy bulge blue optical magnitude, 
and a standard B — R ~ 1.5 to convert to elliptical galaxy 
blue magnitudes to red, one obtains 



i 0fr p FR 

LVJ S J rad,rnod 



-0.65 Mr + 11.1 

+ log(e_ 3 C-iar 2 7/10 "i 2/ 3 5 '7-ii) (19) 



Given the uncertainties in the observed slope, in the scal- 
ing of the coronal density, and in the parameters e, 
and, especially, 77, the good numerical agreement between 
equations ( |l8| ) and ( |l9| ) is rather fortuitous, of course. The 
important points to note are as follows: 

• The basic trend that brighter galaxies, with their 
more massive black holes, should undergo the mag- 
netic switch at a higher jet power is a general feature 
of the model (see figure ||). The ultimate source of 
this trend is that the critical jet power is proportional 
to the size of the ergosphere squared times the disk 
coronal density there. If this density were uniform 
for all black holes, it would produce a trend slightly 
steeper than P™ oc L 2 pt , because the galaxy M/L 
ratio increases with galaxy L opt ( |Faber et al 1987| ); 
however, as densities tend to decrease as the size of 
the accreting system increases, the trend should be 
less steep than this. 

• Although the uncertainties are great, the intercept is 
consistent with the values of the parameters chosen 
here, in particular with the very low value of 77. That 
is, the observations are entirely consistent with the 
jet material originating in the coronae of the black 
hole accretion disks, and grossly inconsistent — by 
more than ten orders of magnitude — with the jet 
ejecting a large fraction of the accretion disk itself. 



Incidentally, the present model is not inconsistent with 
observ ations that F R I jets may decelerate on kiloparsec 
scales (Laing 1996). Indeed, it suggests a reason for such 
deceleration to occur more readily in FR I sources. Jets 
with Lmhd < L cr it are transonic in the simulations here 
and in MEGPL. Such jets ar e more prone to deceleration 
than highly supersonic flows ( Bickncll 1985 ). Note that de- 
celeration of relativistic parsec-scale jets to nonrelativistic 
flow on kiloparsec scales is needed in any model, magneti- 
cally switched or not, to explain head-tail and wide- angle- 
tail sources in clusters, whose jets are bent by motion 
of only a few thousand kms -1 through the intracluster 
medium. 

4.2. Black Hole Spindown and the Early Evolution of FR 
II Radio Sources 

Because the black hole spin is the assumed source of 
power, there should be a critical spin j cr a — correspond- 
ing to a critical hole angular velocity Qn,crit — such that 
the magnetic switch is triggered when the hole is spun up 
to this value. Equating P^ a d, mo d and Prad,mod (or L MHD 
and L cr it) one obtains 



Jcrit 



n ori -2/5 1/10 . -1/5 1/2 

0.20 a_2 rn 9 w_3 V-11 



(20) 



Now, while this relation is a very weak function of mg and 
m, and a is considered to be known to within a factor 
of ten or better, the exact value of j cr a is still very un- 
certain because of the parameter rj. Indeed, a large value 
(dense corona) would make j cr n > 1, implying that even 
a maximal Kerr hole would not trigger the switch. Below 
observational tests are proposed to confirm the low choice 
of rj in this paper. 

Note that, independent of the magnetic switch model, 
if the source of the radio jet power is, indeed, the black 
hole rotational energy, then simple arguments show that 
at the FR I/II transition, j must be of order 0.1 in any 
case. From the OLWL results and the bivariate radio lu- 
minosity function itself, it appears that, for a galaxy of a 
given mass, the maximum possible radio power is about 
100 times the power at the FR break. Therefore, for any 
model where P ra d oc j 2 , we find that 



Jobs 



pFB, 
rad 

prnax 
rad 



1/2 



0.1 



(21) 



The magnetically-switched, rotating black hole model 
further predicts that FR II sources contain the seeds of 
their own destruction. The jet extracts rotational energy 
from the black hole until j < j cr it, whereupon the jet tran- 
sitions to an FR I source. Now, in the absence of external 
effects, the spindown rate is proportional to j , and there- 
fore the hole rotational energy 



E rnt = Mc z 



(1 + J 2 /4) 



1/2 _ 



1.6 x 10 62 ergm 9 f 



(22) 



rendering the spindown an exponential process, with a 
time scale 



E r 



Tspindown 



4.1 x 10 8 yr oT 2 1/10 1/10 ml 3 4/5 Q 1 (23) 
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where £ is the duty cycle parameter defined earlier. This 
spindown time is ^ 10 9 yr time scale found for the cosmic 
evol ution of powerful cxtrag alactic radio sources ( [Schmidt 



1972 fc|Wall fc Jackson 199 



1 



A plausible scenario, then, 
is that a fraction of black holes are born with very high 
spin (FR II sources) and subsequently experience a strong 
density evolution, disappearing after T sp i n d OW n or so. Note 
that the local FR I space density exceeds that of FR II 
sources at all epochs, so that expired FR II sources could 
easily be hiding in the FR I population ( Jackson fc WaD 
1999|Jjackson, private communication). 
Note also that, like the accretion time 



M I M 



4.5 x 10 12 yr mZ\ Ci 



(24) 



Tspindown is virtually independent of black hole mass (al- 
though much shorter). Therefore, in the absence of other 
effects, even stellar mass black holes like GRO J1655-40 
or GRS 1915+105 are expected to retain their spin over 
~ 10 9 yr as they use their rotational energy to drive radio 
jets. Such a long spin lifetime should have implications for 



mode )s for galactic microquasars. 



During the entire lifetime of the radio source, an en- 
ergy comparable to E rot should be liberated into the lobes. 
That is, extended radio sources are the remnants of for- 
merly rapidly-rotating black holes. This amount of energy 
is consistent with, and in fact somewhat larger than, the 
total energy content of extended ra dio sources in r elativis- 
tic particles and magnetic fields ( DcYoung 1976| ). This 
allows for a significant fraction of this energy to have been 
converted to other forms {e.g., heating of the intergalactic 
medium) or radiated away. 

4.3. Little Evolution for FR I (and FR II) Sources at 
Late Times? 

The fact that E rot is proportional to j 2 , and not to just 
j, for example, means that the spindown time is the same 
no matter how rapidly or slowly the hole is rotating. The 
model predicts therefore, that after turning into an FR I 
source, the object will continue to display pure luminosity 
evolution on a T sp i n down time scale. However, rece nt anal- 
yses of the radio source counts and rcd shift data (Urry & 



which, with efficiency factors similar to those used 
above, could generate radio emission of up to ~ 
10 21 WHz -1 . A radio source of this power is not 
expected to evolve much, unless m evolves. How- 
ever, this is a very weak power — much smaller than 
many FR I radio galaxies, which also presumably 
do not evolve much. Furthermore, if the accretion 
events have rather random angular momenta, then 
this minimum jet has no preferred direction. It is 
likely to point in random directions with each event, 
never developing into an extended source, and per- 
haps appearing radio quiet. 

• Periodic spinup of the black hole. On the other 
hand, if the short-term accretion events have sys- 
tematically the same angular momentum direction, 
as might be the case for a galaxy collision or merger, 
then while they would keep the above weak jet ori- 
ented in the same direction, the y also would tend 
to re - orient and spin up the hol e (Wilson & Colbert 
199?;; Natarajan fc Pringle 199q ). This, or the direct 



merger of two black holes, would re-kindle the ergo- 
spheric jet. Then, if the galaxy merger rate evolved 
on a very long time scale (say, ^ 5 x 10 9 yr), the 
shorter time scale spindown process would become 
slave to the spinup process, with the average radio 
power being a function of the rate at which angular 
momentum is accreted onto the hole. At low red- 
shifts, then, both FR I and the occasional re-kindled 
FR II objects would evolve on the longer merger evo- 
lution time, not T spin down- 

4.4. Radio Loud vs. Radio Quiet? 

The model proposed here deals mainly with the differ- 
ence between two radio loud populations — the FR I and 
II sources. However, it also has implications for the dif- 
ference between radio loud and quiet quasars as well. In 
particular, it suggests that radio quiet objects should be 
identified with very slowly-rotating, nearly-Schwarzschild 



black holes. At first, the minimum jet power (equation 
Padovam 1995] |WaU fc Jackson7997D' conclude tEaFEEc H appears to be a problem with this suggestion, as it 



FR I population evolves little or not at all. Modeling this 
sort of behavior is very difficult, especially in a universe 
where everything involved with AGN seems to evolve on 
at least a cosmological time scale (galaxy formation and 
mergers; the rate of stellar death [which may feed gas into 
the nucleus]; even radio quiet quasars; etc.). Below are 
some ideas for halting, or at least slowing, the evolution 
of FR I sources in this model; however, much more work 
needs to be done on this problem: 

• Constant low-powered jet from outside the ergo- 
sphere. After spinning down, the black hole loses 
most of its ergosphere and its ability to create a large 
poloidal field strength in the disk. The field st rength 



predicts that all active galaxies should produce some sort 
of extended radio emission. 4 However, even if an AGN 
does produce an extended radio source power of 



Tjmin 
r ad 



< 9 x 10 20 WHz 



9/10 

V 



£-3' 



-1/10 . 6/5^ 



(26) 



is expected to drop to the value suggested by Livio 

yielding an expected jet 



Ogilvie 
power of 



j mm 
^MHD 



Pringle (1999) 



9 x 10 40 ergs" 



a 



1/10 9/10 . 6/5 



111 



111 



(25) 



the apparent flux would be quite small: less than a few /iJy 
for Seyfert galaxies at z > 0.03 and quasi-stellar objects at 
Z > 1, and a few millijanskies total for nearby spiral galax- 
ies. This amount of extended emission in spirals could be 
hidden easily by emission from supernova remnants, etc.. 
Therefore, there should be a minimum spin rate j m in for 
the hole such that when j < j m i n , it will appear radio 
quiet. While detailed simulations of jet production from 
rotating black holes, plus observational selection factors, 
need to be taken into account to determine a precise value 



4 Of course, as discussed earlier, due to a lack of a consistent long-term jet direction, this power may never manifest itself as an extended 
radio source. Here it will be assumed that it does; if not, the conclusions will be stronger. 
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of jmin, we can estimate it here by equating the MHD 
power generated by the ergospheric and Keplerian disk 
jets (equations E2I and |2j3) and solving: 



Jr, 



= 0.83 x 10 



-3 -1/10 -1/10 • 1/5 

a n TOn TO q 



(27) 



For spins slower than this, the wind/jet from the Keplerian 
portion of the disk outside the ergosphere will dominate 
the output MHD power. 

Even if the minimum jet power turns out not to be a 
problem, one still needs to explain — quantitatively as 
well as qualitatively — why the radio power distribution 
is bimodal: radio quiet objects appear to be a separate 
class, no t merely the fainter versi on of the powerful radio 
sources ( Kellermann et al. 1989 ). One possible solution 
to this problem is the re-kindling the radio source thr ough 



galax y mergers or co llisions. This is not a new idea (Wil- 
son & j Colbert 1995^ but the model presented here, es- 



pecially equation (23), gives a physical and quantitative 
basis for it. In the early universe, holes would have been 
born with a spectrum of rotational speeds. Those rotat- 
ing rapidly would have produced extended radio sources, 
including holes in some (pre-)spiral galactic bulges, while 
those rotating slowly would have been radio quiet. All of 
these objects should have spun down in a few T sp i n down 
times and become radio quiet, according to the model. 
However, a subgroup — those in galaxies that experience 
repeated collisions and mergers — will be spun up again, 
at least partially to FR I class and very occasionally to FR 
II strength, re-kindling the radio source. This scenario is 
consistent with powerful radio sources appearing predomi- 
nately in elliptical galaxies, which now are believed to form 
from, or be strongly influenced by, mergers. An important 
corollary is that, early in their history, even the bulges of 
a significant fraction of spiral galaxies should have under- 
gone a one-time extended radio source phase, lasting a few 
Tspindown times and producing both FR I and II morpholo- 
gies, but with considerably weaker luminosities than those 
associated with present-day giant elliptical radio sources. 

The radio loud/quiet problem thus may be one of envi- 
ronment rather than being intrinsic to the system. How- 
ever, before adopting such a simple model, it will be im- 
portant to understand the strange case of the microquasar 
GRO J1655-40, which can be cit her radio loud or r adio 
quiet during its X-ray outbursts ( Tavani ct al. 1996 ). It 
is not clear yet whether this is simply a short-term phe- 
nomenon or is fundamentally related to the macroquasar 
radio loud/quiet problem. So far, no model, including the 
one in this paper, appears able to explain completely this 
object's behavior. 

4.5. Secondary Evolution in Accretion Rate and the Lack 
0/ FR I Quasars 

It has been assumed in this paper that on a spindown 
time the average to is fairly independent of black hole mass 
and remains at roughly ~ 10 -3 . However, on a cosmologi- 
cal time scal e this is not likely to be the case. Indeed, there 
is evidence ( [Wall fe Jackson 1997| ; [lackson fe Wall 1999| ) 
that some sources (class 'A') have a high accretion rate, in- 
dicated by strong, broad line emission, while others (class 
'B') have a low (perhaps ADAF-like) rate. For the case 
where the latter objects are core-dominated and appear- 
ing as BL Lac objects, Wall & Jackson note that at least 



some must have FR II sources as their parent population, 
not just FR I sources. One possible explanation of this 
is that the average to may have been much higher in the 
early universe when the merger rate was high, and lower, 
perhaps even ADAF-like, now. Indeed, it appears that it 
i s only the class 'A' o bjects that display rapid evolution 
( Jackson fe Wall 1999 ). Then, present-day re-kindled FR 
I and II sources would both have a low average accretion 
rate, as well as a low evolution rate, and both have the po- 
tential for appearing as a BL Lac object if viewed end-on. 
Not all FR II sources seen this way would have to appear 
as quasars. 

Ascribing a high accretion rate to class 'A' objects solves 
another puzzle concerning quasars — the fact that few, if 
any, radio loud quasars display FR I morphology. For 
to ~ 1, equation (j2Tj) gives j cr n ~ 0.05, implying that 
even slowly-rotating black holes would produce an FR II 
morphology when the accretion rate is high. The right 
panel of figure |^ illustrates that the FR I portion of the 
quasar diagram is significantly smaller than that for radio 
galaxies; the exact size of the region depends critically on 
the definition of radio quietness. 

5. TESTS OF THE MODEL 

5.1. Observational Tests 

This paper has shown that the magnetic switch mech- 
anism, triggered by a rapidly spinning black hole, is con- 
sistent with a variety of observations concerning the FR 
I /II break. However, Bicknell (1995) has shown that an 
ISM model is consistent as well. In order to distinguish 
between these two, the following observational tests are 
proposed. The current model predicts the following: 

1. The ensemble average of VLBI jet speeds for sources 
with FR I morphology should be distinctly slower 
than that for FR II sources. Until detailed general 
relativistic MHD simulations of very tenuous coro- 
nae in Kerr black hole ergospheres can be done (c/. 
Koide et al. 1998b| ), the exact values of the two mean 
velocities will be uncertain. The FR II mean should 
be highly relativistic (i.e., 7 ^ 5 — 10 or greater; 
see MEGPL). However, the FR I value may be more 
than mildly relativistic itself, as the escape velocity 
from the ergosphere is very close to c. Note the em- 
phasis on FR morphology, rather than radio power. 
Radio power should not be used to distinguish be- 
tween FR class in this test, unless it can be shown 
(using the OLWL relation) that the galaxy in ques- 
tion lies well within either the FR I or FR II region 
of the radio-optical diagram. Faint galaxies can have 
low-power radio sources that still have FR II mor- 
phology, and vice-versa. 

2. Direct measurement of the accretion disk corona 
density in the ergosphere, using X-ray observations 
for example, should yield rather low mass densities 
(~ 10~ 15 gcm~ 3 ), perhaps either a 10 12 cm' 3 e + e~ 
pair plasma or a very tenuous 10 9 cm~ 3 p + e~ corona. 
These low values are needed for L cr n to be such that 
a radio power of P^ mod ~ 10 25 WHz _1 is gener- 
ated at the FR break. Significantly higher measured 
densities would render the corona too heavy for the 
magnetic switch to be triggered by any reasonable 
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jet power, and, in this model, always would predict 
slow jet velocities. 

3. Recently, methods have been developed to estimate 
the spin of black holes by measuring the rotati on rate 
of the innermost part of their accr etion disks (Zhang 
et al. 1997|; |Sobczak et al. 1998|). If such methods 



can be applied to extragalactic radio sources, one 
should find that the ensemble average of black hole 
spins for sources with FR I morphology should be 
distinctly lower than for those with FR II morphol- 
ogy- 

4. There should exist a population of high-redshift, 
extended radio galaxies associated with the bulges 
of spiral or pre-spiral galaxies. They should in- 
clude both FR I and II morphologies, but be weaker 
copies of the giant sources associated with ellipti- 
cals. The radio sources themselves should be read- 
ily detectable (for mg = 10~ 2 and a redshift of 3, 
Prad ~ lO^-^WHz" 1 or S ~ 4 - 4000/^Jy), but 
their optical counterparts would be very faint if un- 
dergoing an initial burst of star formatio n (Mr ~ 25 
for a 10 10 Mq bulge at a similar redshift; |Meier 1976 ) 
and even fainter (Mr ~ 30) if not. Preliminary es- 
timates of their numbers are on the order of a few 
hundred to a few thousand per square degree — con- 
siderably smaller than t he larger number of star burst 
galaxies at faint fluxes (Windhorst ct al. 1995). 

The first and third of these tests would provide support 
for the magnetic switch mechanism for the FR I /II tran- 
sition. The second would support the choice for a small 
rj. The third and fourth provide evidence for the proposal 
that the black hole spin is the switch triggering mechanism 
and the ultimate source of the radio power. 



6. CONCLUSIONS 

The magnetically-switched, rotating black hole model 
for the FR I/II break is consistent with current black hole 
accretion models for active galactic nuclei, with MHD sim- 
ulations of jet production, and with relativistic wind the- 
ory. It predicts that a sharp transition in jet morphol- 
ogy should occur with increasing radio power, as observed, 
and that the high luminosity jets should be fast and ki- 
netic energy dominated. It yields the correct trend of this 
transition with galaxy luminosity (larger galaxies should 
transition at larger radio power), and is consistent with 
the actual value of the FR I/II break radio power, if the 
average density of accretion disk coronae in black hole er- 
gospheres is rather low (~ 10~ 15 gem -3 ). The model sug- 
gests a long-term memory mechanism — the black hole 
spin — for assuring that radio sources maintain the same 
morphology over periods of time that are at least as long 
as the flow time from galactic center to lobe (1 6 ~~ 7 yr). In 



that regard, it builds on an earlier suggestion ( [Rees 1978 ) 
that black hole spin is the mechanism determining the di- 
rection of the radio jet (although the black hole spin itself 
may ultimately be ti ed to the angular momentum of merg- 
ers or collisions; see Natarajan & Pringle 1998). Because 
black hole rotational energy drives the radio jet, the model 
predicts spindown and, consequently, cosmic evolution of 
the FR II space density on a time scale of order 10 9 yr 
or less. It also suggests a possible model for radio quiet 
quasars (Schwarzschild black holes) and for why radio loud 
quasars display predominantly FR II morphology. Finally, 
the total amount of rotational energy liberated by the jets 
from a rapidly-rotating hole is more than adequate to ac- 
count for the entire energy content of a typical extended 
radio source. 

Several observational and numerical experiments have 
been proposed to test this model. The most important of 
these is the VLBI jet velocity test, in which objects with 
FR I morphology should have a distinctly slower ensemble- 
averaged jet speed than that for FR II sources. Other tests 
involve a similar effect for the black hole spin, a predic- 
tion for the typical density of a black hole accretion disk 
corona, and a search at high redshift for a population of 
faint extended radio sources associated with spiral or pre- 
spiral bulges. Proposed numerical experiments include 
general relativistic simulations confirming that the mag- 
netic switch process works and that output MHD power 
depends critically on black hole spin. 

Note added in proof: In discussing the literature on 
the expected power output of a rotating black hole, the au- 
thor omitted an important reference. P. Ghosh & M. A. 
Abramowicz (1997, Mon. Not. Royal Astron. Soc. 292, 
887) also discussed this subject at length, assuming, as is 
Romanova et al. (1998)| , with closed coronal loops, done here, that B p0 w B d0 , and obtained expressions for 

the MHD power output of a rotating hole similar to those 
in the present paper. (Note that these authors took all 
of this power to be due to the Blandford-Znajek process 
occurring very near the horizon.) However, the absolute 
magnitude of their power, is about two orders of magni- 
tude smaller than that in equation (12). Part of this dif- 
ference is due to a factor of 8 in the definition of the MHD 
power (which was discussed by these authors) , but most of 



5.2. Numerical Tests 

A number of numerical tests will be almost as important 
as the observational ones: 

1. Firstly, as suggested by MEGPL, the magnetic 
switch needs to be confirmed by other groups per- 
forming numerical simulations, especially with a 
general relativistic MHD code, even if just for 
Schwarzschild black holes. 

2. Secondly, general relativistic MHD simulations in 
Kerr geometry are needed to confirm that a substan- 
tial jet, whose power is directly related to the black 
hole spin, can be produced simply by threading the 
ergosphere with a magnetic field that is anchored 
further out. A simple vertical field should be suf- 
ficient initially, but simu lations similar to those of 



would provide a more realistic test. 

3. Finally, it will be important to reproduce in Kerr 
geometry the simple simulations performed in this 
paper, showing that increasing the black hole spin 
eventually triggers the magnetic switch. It also will 
be important to determine the speed, Mach num- 
ber, and collimation properties of the outflow when 
the black hole spin is well above fl cr it to see if the 
flow de-collimates (as it did in these non-relativistic 
simulations) or remains well-collimated. 
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it is due to the present paper recognizing that the accretion 
disk can extend well interior to 6GM/c 2 and even inside 
the ergosphere for rapidly rotating black holes. This yields 
much larger magnetic field strengths near the static limit 
and black hole horizon. In fact, one way in which Ghosh 
& Abramowicz (1997), Livio et al. (1999), and this paper 
can be brought into rough agreement is to 1) extend the 
disk into the ergosphere (which would increase the power 
computed by Ghosh & Abramowicz) and 2) treat that disk 
as an advection-dominated one with Hq/Rq ~ 1, yielding 
B p o « Bdo (which would increase the power calculated by 



Livio, Ogilvie, & Pringle). 
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APPENDIX 



METRIC-SHEAR-DRIVEN MHD DYNAMOS 

This appendix shows that shear in a general relativistic metric, like that displayed in Kerr geometry, can serve to 
enhance the local magnetic field strength. Magnetized plasma accreting onto a rotating black hole, therefore, will extract 
energy from the hole, converting it into enhanced local magnetic field strength. The effect occurs even in absence of shear 
in the plasma fluid flow itself and even if the magnetic field lines do not thread the black hole. 

In a classical MHD flow the magnetic field evolution equation, derived from Faraday's and Ohm's laws with infinite 
conductivity, is 



dB 

~dt 



v-VB = B-Vv BW-v 



(Al) 



where B and v are the magnetic field and velocity vectors. For an initially poloidal field and rotational flow only, the 
component of this equation becomes, in spherical-polar coordinates, 



dt 



r sinO 



B r 



on 

dr 



B„ 



i on 

r ~d9 



(A2) 



where = v$/(r sinQ) is the angular velocity field of the flow. Thus, if there is shear in a classical rotating magnetized 
fluid flow, a poloidal magnetic field will serve as a seed for generation of a potentially large azimuthal magnetic field 
that gro ws with time. T his "omega-effect" is a necessary (th ough not suffic ient) condition for an "alpha-omega" dynamo 
to exist (Roberts 199S). According to Cowling's theorem (Cowling 1934), for a full dynamo, there also must be non- 
axisymmetric poloidal fluid motions (due to turbulence, for example) that enhance the poloidal magnetic field as well. 
Such turbulent motions arc thought to be a com mon occurrence when the rotating shear flow occurs in a gravitational 
field (3tone et al. 1996; Brandenburg ct al. 1996). So, the presence of an omega-effect in a gravitational field is a strong 
indication that an alpha-omega dynamo may be operating and generating a substantial magnctosphcrc 



The theory of gravito-magnetic dynamos in the Kerr metric has been developed by Khanna & Camenzind (1996a) 
raising the possibility that even an axisymmetric flow can have growing, sclf-cxcitcd modes. Although such modes so 



far have not been found in kinematic numerical simulations (Brandenburg 1996; Khanna fc Camenzind 1996b ), the non- 
validity of Cowling's theorem in a general relativistic metric has been analytically confirmed by [Nunez (1997) , who showed 
that there are growing axisymmetric modes when steep gradients exist in the plasma angular velocity £1 relative to the 
rotating space. Here, however, we consider the simpler situation where the plasma flow itself does not shear, but there still 
is shear due to the general relativistic effect of frame-dragging by a rotating black hole. While this does not demonstrate 
the existence of a dynamo in the system, it does satisfy several of the necessary conditions for such — in particular, an 
omega-effect occurring in a gravitational field. 

To compute the evolution of the azimuthal magnetic field, we us e two different approaches — one that uses the covariant 



four-dimensional formalism of Misncr, Thornc, fc Wheeler (1973) (MTW) and one that uses the "3+1" decomposition of 
Thornc, Price, fc MacDonald (1986) (T PM). Both give ins ight into the magnetic field enhancement process. The general 
relativistic analogy of equation (Al) is ( Lichncrowicz 1967) 



wB -g — B H u - } p — B u -3 



(A3) 



which uses the standard MTW conventions (Einstein summation, semicolon to signify covariant differentiation, and comma 
for ordinary differentiation) . Expanding the covariant derivatives and solving for evolution of B a 



/ B a i 0=B^u a ,3 + s^(r%3 - TVr) - 
= B^u a t8 + BPvFcfo" - B a u 3 .3 



B a vP-i- 



(A4) 



where the r™^ 7 are the connection coefficients and the antisymmetric commutation coefficients are given by the components 
of basis vector commutators 

C /3 7 Q = [ e p, e iT (A5) 
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The metric shear C/3 7 a vanishes for B a expressed in a coordinate basis. But it does not for a physical, locally-Lorentz, 
orthonormal frame (one in which local observers measure physical quantities, and in which magnetic field has the units of 
Gau ss for all components). In such an orthonormal frame, where we denote the coordinates with a "hat" accent, equation 
(A5) becomes 

where D 1 a is a coordinate transformation that diagonalizes the metric, converting to a locally-Lorentz system, with inverse 
CP 

To compute Cg~ Q , and therefore the evolution of the magnetic field for rotating black hole magnetospheres, we use the 
zero-charge Boyer-Lindquist version of the Kerr metric (MTW) written in the following form 



ds 2 = 



■di 2 



^dr 2 
A 



p 2 de 2 



<dty 



(A7) 



where the lapse function and angular velocity of frame dragging are given by 

2aMr 



a = — A 1 / 2 
E 



the various geometric functions are given by 

A = r 2 - 2Mr + a 2 
S 2 = (r 2 +a 2 ) 2 - a 2 Asin 2 9 



p — r + a cos 
S 



sinQ 



P 



a — J/M = jM is the specific angular momentum, and physical quantities are now in geometric units (c = G = 1). 
This metric can be diagonalized with many transformations, of course. We choose one that rotates with the same angular 



veloci ty u as the frame-dragg ed space. This is the so-called Zero Angular Momentum Observer (ZAMO) system of Thornc 
Price, fc MacDonald (1986)|. The non-zero elements of this transformation and its inverse are 



*~ t 


= l/a 


£ r " 


= A 1 / 2 / p 




= uj 1 a 


ch 


= a 




= p 1 A 1 / 2 




= — WU) 



£ 6 e = 1/P 



re 



P 
W 



with the ef and eg axes parallel to e r and eg, respectively. 

To illustrate that the metric shear acts like a fluid shear flow, consider the situation where there is no plasma motion 
with respect to the orthonormal system, so that the fluid four-velocity is given by u a = (— 1, 0, 0, 0). Expressing equation 



(A4) in the orthonormal ZAMO frame, the evolution of the (f> component of the magnetic field, using equation (A6), is 



B" 



A 1 / 2 

■w 

a p 



or 



p § 1 9uj 
A 1 / 2 89 



C (j>i 



(A8) 



with the time and azimuthal parts vanishing. Comparing with equation ( A2 ) , then, we see that there is an exact analogy 
with the classical rotating fluid shear case. Any poloidal magnetic field line penetrating the space near the rotating black 
hole is sheared into azimuthal field at the expense of the rotational energy of the black hole, whether or not the magnetic 
field line actually threads the hole. This is independent of any shear in the fluid flow itself that may or ma y not be present. 

Equation ( [Asf ) also can be derived using the "3+1" language of [Thornc, Price, fc MacDonald (1986] (TPM). In that 
case the metric equation (A7) is expressed in terms of the lapse and shift functions and a 3-metric 



ds 2 



-a 2 dt 2 + -fij(dx l + I3 l dt)(dx j + (3 j dt) 



(A9) 



where 7y raises or lowers indices on the shift 3-vector (3 l and 3-velocity v % . Comparing (A7) and (A9) gives /3 l = (0, 0, —lo), 
and components of (the diagonal) 7^ can be determined similarly. Physical quantities like magnetic field and velocity 
vectors are measured in a reference frame that rotates with the ZAMOs, but they and the gradient operator V are 
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expressed in terms of absolut e (B oyer-Lind quist) spatial and time coordin ates. Under these conditions, the general 
relativistic analog to equation (Al) is (TPM; Khanna & Camenzind (1996a)) 



dB 

~dt 



- £ P B + av-VB = B-V(av) - BV-{av) 



(A10) 



where £p represents the Lie derivative along /3 

£pB = (3 VB B V/3 

and the frozen-held condition E — —vxB and solenoidal condition V ■ B = 
situation where there is no fluid motion with respect to the rotating space 



(All) 

: have been applied. We again consider the 
v = 0). Furthermore, because the gradients 



now are expressed in a coordinate frame, the shear terms in equation (All) vanis h, an d the Lie derivative becomes 
simply £pB l = ft B l _j — B J ft ^ As a result, the azimuthal component of equation (A10) becomes, in Boyer-Lindquist 
(un- hatted) coordinates, 

B*, t - P*B*4 = B j u^ (A12) 
where j = (r,9) only, since d(3/d(f) = —duj/d<p = 0. The left-hand side is proportional to the time derivative in the 



ZAMO frame (= aB^ f/w), so transformation of equation ( A12) to the locally-Lorentz orthonormal coordinates yields an 



equation identical to (A8). That is, even though the 3+1 formalism is expressed in the Boyer-Lindquist coordinate frame, 
the metric shear still is present; but that shear is now embodied in the shift vector f3, rather than in the commutators 

C/3 7 ". 

Since the sheared space near a Kerr black hole is so similar in nature to that of a sheared accretion disk (i.e., rotating 
shear flow in a gravitational field), it is likely that an instability similar to the Balbus-Hawley one also is operating. Thus, 
there should be not only an omega-effect enhancing the azimuthal magnetic field, but also an alpha-effect, due to the 
magneto-rotational-generated turbulence, enhancing the poloidal field as well. Since this process draws its energy from 
the hole's rotation, the strength of the magnetosphere is expected to continue to grow long after all the accretion disk's 
angular momentum has been removed, and to levels far greater than could be supported by the accreting matter alone. 
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Fig. 1. — Results from 13 MHD simulations that are similar to those in MEGPL, but with angular velocity of the injected material used 
as the independent parameter rather than poloidal magnetic field strength. Left panel: jet velocity (filled squares) and Mach number (open 
circles) vs. angular velocity, normalized to the Keplerian value. Right panel: power output of the disk system in the form of jet kinetic power 
(filled squares), total kinetic power (wind plus jet, open diamonds), and total power in all forms (filled circles). For small values of f2/f2x, 
the disk power is dominated by advected magnetic, thermal, and kinetic energy flowipg into the corona (whose sum varies little or not at all 
with f2). For large values of Q/flpc, the disk's magneto-rotational power (equation W) dominates. The predicted value of Q, cr it (equation ph 
has been modified to account for the dynamical adjustment in conditions that occurs in the injection region when Q ^ Q.^ (see text). When 
Q < Cl cr it, the total power is dominated by advected magnetic energy, and the kinetic energy is dominated by that in the looscly-collimated 
wind. As Q exceeds fl cr n, the velocity and Mach number increase by an order of magnitude or more and the jet kinetic power becomes a 
substantial fraction of the total. In these non-relativistic simulations, eventually rotation becomes so fast that the flow de-collimates into a 
mostly radial wind. 



Fig. 2. — The theoretical radio-optical plane for sub-Eddington accreting radio galaxies (left panel, m = 10 — 3 ) and near-Eddington 
accreting quasars (right panel, rh = 1.0). For galaxies, the optical luminosity scales roughly with the galaxy mass, as does the black hole 
mass. For quasars, the optical luminosity scales linearly with the black hole mass. In both cases, the radio power scales with the square of 
the normalized black hole spin j. Note that only the extended (unbeamed) radio power is considered here; orientation effects are ignored in 
these diagrams. The upper solid lines show equation (h4) for a maximal Kerr hole. Middle solid lines show thi^theoretical FR I/FR II break 
(equation 171), which agrees well withJiie observed OLWL relation for radio galaxies (c/. equations ha and |l9| ). Lower solid lines show one 
measure ofradio quietness (equation Brjh, below which the Keplerian disk MHD power dominates the ergospneric disk MHD power. Note 
the much smaller region spanned by FR I quasars than by FR I galaxies. FR I sources should comprise a much smaller fraction of quasars 
than of radio galaxies. The right panel also shows another estimate of radio quietness that makes the FR I quasar region even narrower — 
dashed lines where the radio to optical/UV flux ratio (1 GHz to 1 PHz) attains values of 10 and 100. The optical/UV flux was estimated by 
integrating the power output of a standard Shakura & Sunyaev accretion disk outward from the radius where the disk achieves an effective 
temperature of 10 5 K. The vertical dot-dashed line shows roughly where the optical/UV luminosity is equivalent to the stellar luminosity of a 
typical spiral galaxy, dividing Seyfert galaxies from quasars. As discussed in the text, the great majority of radio loud Seyferts are expected 
to occur only at high redshifts. 
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